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Experimental studies of power corrections with e + e~ data are reviewed. An overview of the available data for jet and 
event shape observables is given and recent analyses based on the Dokshitzer-Marchesini- Webber (DMW) model of 
power corrections are summarised. The studies involve both distributions of the observables and their mean values. 
The agreement between perturbative QCD combined with DMW power corrections and the data is generally good, 
and the few exceptions are discussed. The use of low energy data sets highlights deficiencies in the existing calculations 
for some observables. A study of the finiteness of the physical strong coupling at low energies using hadronic r decays 
is shown. 

\o 
o 
o 

<N 1. INTRODUCTION 

\ Data from hadron production in c + c annihilation experiments has been one of the main sources of experimental 
information for the study of non-perturbative QCD of hadron jets, see e.g. 0,0,0,0- Compared to other processes 
where power corrections have been studied such as hadron production in deep inelastic lepton-nucleon scattering |2| 
there are two main advantages: 

> ; 

Leptonic initial state The incoming electron and positron are leptons and thus there is no interference between 
initial state and the hadronic final state. This makes the interpretation of of the final state in terms of QCD 
■ much more robust. 

Large range of energies The currently available data from e + e~ annihilation to hadrons span a range of centre- 
of-mass (cms) energies of more than an order of magnitude. Since many effects in QCD depend on the energy 
^ ' scale comparing data sets recorded at different cms energies yields strong tests of the theory. 

From reconstructed hadronic final states jet production or event shape observables are calculated, see e.q. |1|, 
As an example the observable Thrust T is defined by T = max^(^\ |n • \Pi\) where pi is the 3-momentum 

of particle i and \ri\ — 1. A value of T = 1 corresponds to an ideal 2-jet event, events with three final state objects 
have 1 > T > 2/3 while completely spherical events with many final state particles can approach T = 0.5. Event 
shape observables y are conventionally defined to have y — for ideal two-jet configurations and thus the Thrust is 
often analysed in terms of 1 — T. The observables can be classified as 3-jet or 4-jet observables depending on whether 
y > occurs in the cms for final states with > 3 or > 4 objects. 



1.1. Experiments 

The energy of the electron and electron beams in most e + e~ annihilation collider experiments is equal due to the 
simplified construction and operation of a circular accelerator with only one beam line. As a consequence the cms 
of hadronic final states produced in e + e _ annihilation coincides with the laboratory system in the absence of hard 
initial state radiation (ISR). 

A typical e + e~ annihilation experiment is shown in figure^ The JADE detector (for details see 5]) was operated 
at the PETRA accelerator Q from f 979 to f 986 at cms energies 12 < ^/s < 46.8 GeV. The design of the experiment 
is similar to the more recent LEP experiments. The interaction region, where the electron and positron beams collide, 
is surrounded by gaseous tracking detectors (Jet Chamber) inside a solenoidal magnetic field of 0.48 T. The tracking 
detectors are inside a calorimeter consisting of 2712 lead glass blocks to measure the energy of charged and neutral 
particles. Tracking and calorimeter efficiency is good in the region | cos#| < 0.97 were 9 is the angle w.r.t. the beam 
axis. Therefore most e + e~ annihilation events arc fully contained in the detector and complete reconstruction with 
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high efficiency and good precision is possible. This reduces the size of experimental corrections for acceptance and 
resolution and thus the corresponding experimental uncertainties. 
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Figure 1: Drawing of the JADE experiment. The electron and positron beam come from the left and right and collide in the 
centre of the experiment. 



1 .2. Overview of data 

The experimental conditions vary significantly between the various cms energies of the data sets. Figure El (left) |7| 
shows the measured cross sections of several e + e _ annihilation processes as a function of the cms energy y/s. The 
\/s dependence as expected from QED at low y/s is clearly visible. At the Z° peak around y/s = 91 GeV the cross 
sections increase by a factor of about 100 compared to their lowest values at y/s ~ 60 GeV or y/s > 130 GeV. 

At y/s < m-zfi background from two-photon interactions with hadronic final states (e + e~ — ► 77), from production 
of r-lepton pairs decaying to hadrons and from hadronic final states with ISR and thus reduced cms is significant. 
Requiring more than four tracks of charged particles (events with exactly four tracks not in a 1- vs 3-prong config- 
uration may also be accepted reduces the backgrounds from r-pairs and two-photon interactions to negligible 
levels. The data collected on or near the Z° peak are essentially free of backgrounds due to the high energy and 
large cross section. ISR effects are also suppressed. The high energy samples taken by the LEP experiments during 
the LEP 2 phase have again backgrounds from two-photon interactions, ISR and for y/s > 2m\y ~ 161 GeV from 
production of four-fermion final states dominated by W + W~ pairs decaying to hadrons. 

As a representative example figure|5] (right) shows the distribution of the event shape observable Thrust T measured 
by L3 at y/s ~ 189 GeV |9(. The simulation of signal and background samples shown in the figure agrees well with 
the data. The selection efficiency is e w 90% and the purity of the s amp les is about p 1=3 80% 0. The remaining 
fraction of background from four-fermion events in the OPAL analysis [10( is 2.1% at 161 GeV and 6.2% at 207 GeV. 

The experimental corrections are typically 0(10)% except at the edges of phase space. They reach values of 50% 
only for the LEP 2 samples because of the reduced selection efficiency of multi-jet final states due to the additional 
cuts to suppress four-fermion events. 

Some 3-jet observables have had QCD predictions in next-to-leading order (NLO) combined with next-to-leading- 
log-approximation (NLLA) predictions (1 — T, Mjj, P>t, Bw, C, y® 3 )- Most power correction calculations are only 
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Figure 2: (left) Cross sections of several e + e~ annihilation processes as indicated Q. The lines represent theoretical predictions. 
The accelerators where the data where recorded are also shown. The variable s' refers to reduced cms energy after ISR. (right) 
Distribution of event shape observable Thrust T measured at ^fs ~ 189 GeV . The lines and hatched regions represent the 
predictions from Monte Carlo simulations of signal and background. 

available for these observables or a subset thereof. The final LEP results are 0, IToL fill while many older 
measurements are summarised in |l3|. Moments up to the order 5 of these and several other observables are also 
available The 4-jet observables D, B^, Ml, T m i n . and were measured by OPAL and partially by L3 



2. TRADITIONAL DMW STUDIES 

The power correction model of Dokshitzer, Marchesini and Webber (DMW) extracts the structure of power cor- 
rection terms from analysis of infrared renormalon contributions jl4j . Earlier versions of the model are \luL . The 
model assumes that a non-perturbative strong coupling ctg p exists around and below the Landau pole and that the 
quantity cto(fii) — J a^ F (k)dk can be defined. The value of /ii is chosen to be safely within the perturbative 
region, usually = 2 GeV. A study of the branching ratio of hadronic to lcptonic r lepton decays as a function of 
the invariant mass of the hadronic final state supports the assumption that the physical strong coupling is finite and 
thus integrable at low energy scales see below. 

The main result for the effects of power corrections on distributions F(y ) of the event shape observables 1 — T, 
Mh and C is that the perturbative prediction Fpt (y) is shifted [la, [lJl |20j : 

F(y) = F PT (y - CyP) (l) 

where c y is an observable dependent constant and P ~ M/ii/Q(ao(/ii) — as) is universal, i.e. independent of the 
observable 0]. The factor P contains the 1/Q = l/\/s scaling and the so-called Milan- factor M — 1.49 for 
ii f = 3 which takes two-loop effects into account. The non-perturbative parameter ao is explicitly matched with the 
perturbative strong coupling as ■ For the event shape observables P>t and B>w the predictions are more involved and 
the shape of the pQCD prediction is modified in addition to the shift [^ij. For mean values of 1 — T, Mh and C the 
prediction is: 

(y) = (y) PT + c y p ( 2 ) 
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For (Bt) and (B\y) the predictions are also more involved due to the modification of the shape of the distributions. 

Figure |3| (left) shows the results of fits to data for mean values of event shape observables measured over a large 
range of cms energies [llj . The perturbative QCD prediction in NLO (y) PT is combined with the power correction 
calculations as explained above. The combined prediction describes the available data well within the uncertainties 
even at low cms energies where non-perturbative effects are large. The numerical results are shown in table [H The 
dashed lines show as a comparison the results of fits of the same NLO QCD predictions with Monte Carlo based 
hadronisation corrections. The data are also well described but the difference in values of cvs^z ) f° r the same 
observables is about 6% fll|| . 
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Figure 3: (left) Mean values of event shape observables as indicated fll]| . The solid lines are fits of NLO QCD predictions 
combined with power corrections, the dashed lines present fits of the same NLO QCD predictions with Monte Carlo based 
hadronisation corrections, (right) Measurements of EEC at various y/s points as indicated |2l|. The solid and dotted lines 
show the fitted QCD prediction including power corrections using only data points within the solid lines for the fit. The 
dashed lines present the same fit after subtraction of the power correction. 



Figure |3 (right) shows data of the event shape observable EEC measured by DELPHI 21] over a large range of cms 
energies. Superimposed are fits of the perturbative QCD prediction in NLO+NLLA combined with power correction 
calculations [22j . The distribution of EEC is sensitive to a higher order non-perturbative parameter referred to 
as a±. A fit with the three parameters as(mzo), «o(2 GeV) and a\(2 GeV) results in ag(m z o) = 0.117 ± 0.002, 
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a (2 GeV) = 0.48 ± 0.05 and a x (2 GeV) = 0.00 ± 0.04. The results for a s (m z o) and a {2 GeV) are consistent but 
the result for «i(2 GeV) is not in agreement with an expectation a.\(2 GeV) = 0.45 from 22] . 





Exp. 


ALEPH [11] 


DELPHI [21] 


L3 [9J 


MPI [23J 




range [GeV] 


12(35)-206 


12-202 


41-206 


12(35)-189 




Q S (m z o) 


0.1207 ±0.0019 


0.1241 ± 0.0034 


0.1164 ± 0.0060 


0.1217 ± 0.0060 


1-T 


a (2 GeV) 


0.539 ±0.011 


0.491 ±0.018 


0.518 ±0.059 


0.528 ± 0.064 




X 2 /d.o.f. 


69/43 


27/41 


18/14 


50/41 




as(m z o) 


0.1161 ±0.0018 


0.1197 ± 0.0122 


0.1051 ± 0.0051 


0.1165 ±0.0043 


M H , P <*o(2 GeV) 


0.627 ±0.020 


0.339 ± 0.263 


0.421 ± 0.037 


0.663 ± 0.095 




X 2 /d.o.f. 


50/40 


10/15 


13/14 


24/35 




a s (m z o) 


0.1148 ±0.0025 


0.1174 ± 0.0029 


0.1163 ±0.0042 


0.1205 ±0.0049 


Bt 


a (2 GeV) 


0.492 ±0.020 


0.463 ± 0.033 


0.449 ± 0.054 


0.445 ± 0.054 




X 2 /d.o.f. 


7/18 


9/23 


9/14 


24/28 




a s (m z o) 


0.1179 ±0.0028 


0.1167 ±0.0019 


0.1169 ±0.0042 


0.1178 ± 0.0025 


B w 


ao(2 GeV) 


0.467 ±0.037 


0.438 ± 0.049 


0.342 ± 0.079 


0.425 ± 0.097 




X 2 /d.o.f. 


11/18 


10/23 


14/14 


10/29 




as(?n z o) 


0.1228 ± 0.0027 0.1222 ± 0.0036 


0.1164 ± 0.0047 0.1218 ± 0.0059 


C 


a (2 GeV) 


0.461 ±0.016 


0.444 ± 0.022 


0.457 ± 0.040 


0.461 ± 0.048 




X 2 /d.o.f. 


17/18 


12/23 


12/14 


18/26 




a s (m z o) 






0.1046 ±0.0124 




D 


a (2 GeV) 






0.682 ± 0.096 






X 2 /d.o.f. 






24/14 





Table I: Results for as(m z o) and Qo(2 GeV) from the analysis of mean values in the DMW power correction model. The errors 
consist of statistical and experimental uncertainties only for ALEPH and DELPHI and include also theoretical uncertainties 
for L3 and MPI. 

Tables HI and ITU summarise all recently published results from studies of DMW power correction using mean values 
and distributions of event shape observables. The values for as^z ) an d olq{2 GeV) are generally consistent with 
each other except for the following observations: 

• The values for <W2 GeV) from Mh or p are comparatively large. This is partially due to the effects of neglected 
hadron masses |24|. The results from DELPHI take hadron masses into account according to |3] and their 
values for cvo(2 GeV) are lower. 

• The results for a>o(2 GeV) from distributions of Bt and Bw from DELPHI appear inconsistent with the other 
results as noted already in |ll| . 

• The results for cvo(2 GeV) are consistent at the level of 20 to 30% but the errors including theoretical effects 
are partially smaller and are about 10% in the case of the MPI analysis of mean values. The Milan factor M 
has a theoretical uncertainty of about 20% jlj(. However, it is not clear if this uncertainty should be reflected 
in differing values of ao(2 GeV) from different observables. 

The correlations between the two parameters as(mzo) and evo(2 GeV) are about —90% from the fits and between 
about —40% and 0% when systematic effects are considered. The results are also shown in figure The negative 
correlation between as(m z o) and ao(2 GeV) is visible in the figure. 

The values of as{mzo) are consistent with a recent determination of as(mzo) from e + e~ annihilation data for event 
shape and jet observables 0: as( TO z°) = 0.1202 ± 0.0010(exp.) ± 0.0055(theo.). The comparatively small statistical 
and hadronisation uncertainties have been included in the experimental and theoretical errors, respectively. This 
result is based on distributions corrected for hadronisation with Monte Carlo models and can most directly be 
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Exp. 


ALEPH [JJJ 


DELPHI [21] 


MPI [23] 




range [GeV] 


91-206 


45-202 


14(35)-189 




as(m z o) 


0.1192 ± 0.0059 


0.1154 ±0.0017 


0.1173 ±0.0057 


1-T 


qo(2 GeV) 


0.452 ± 0.068 


0.543 ±0.014 


0.492 ± 0.077 




X 2 /d.o.f. 


73/47 


291/180 


172/263 




as(m z o) 


0.1068 ±0.0051 


0.1056 ±0.0007 


0.1105 ±0.0040 


M B , P Qo(2 GeV) 


0.808 ±0.185 


0.692 ± 0.012 


0.831 ±0.149 




\ 2 /d.o.i. 


124/42 


120/90 


137/161 




as(ra z o) 


0.1175 ±0.0074 


0.1139 ±0.0016 


0.1114 ±0.0063 


Bt 


q (2 GeV) 


0.667 ±0.137 


0.465 ±0.014 


0.655 ±0.120 




X 2 /d.o.f. 


181/59 


88/75 


92/159 




as(m z o) 


0.1043 ±0.0048 


0.1009 ±0.0016 


0.0982 ± 0.0073 


B w 


ao(2 GeV) 


0.812 ±0.196 


0.571 ±0.031 


0.787 ±0.153 




X 2 /d.o.f. 


76/47 


106/90 


96/132 




as(m z o) 


0.1159 ±0.0062 


0.1097 ±0.0032 


0.1133 ±0.0050 


C 


q (2 GeV) 


0.443 ± 0.056 


0.502 ± 0.047 


0.507 ±0.082 




X 2 /d-o.f. 


83/54 


191/180 


150/208 




as(m z o) 




0.1171 ±0.0018 




EEC 


a (2 GeV) 




0.483 ± 0.041 






X 2 /d.o.f. 




53/90 





Table II: Results for as{m z o) and ao(2 GeV) from the analysis of distributions in the DMW power correction model. The errors 
consist of statistical and experimental uncertainties only for ALEPH and DELPHI and include also theoretical uncertainties 
for MPI. 



compared with the results from distributions shown in table ITTI where all values of as(mzo) are lower. The ALEPH 
analysis 01 makes a direct comparison of results for as{m z o) using DMW power corrections and Monte Carlo 
hadronisation corrections and finds a difference of about 9%. This difference is larger than the commonly quoted 
total error of about 5 to 6% on determinations of as("iz°) from event shape or jet observables. 



3. EXTENDED DMW FITS 

The MPI analysis of distributions of Bw ^ finds large values for x 2 /d.o.f. when y/s < m z o indicating a beginning 
failure of the available calculations to describe the low energy data. For distributions of a similar analysis |25| 
finds large values of x 2 /d.o.f. when data at low y/s are included. 

Figure (left) shows results of simultaneous fits of as(wzo) and «o(2 GeV) to data for distributions of Bw at 
14 < y/s < 91 GeV The solid lines indicate results of extended fits where an additional term An \u(Q)/Q, 

Q = y/s, with a new free parameter An has been added to the DMW prediction. The fit describes the data well 
and yields x 2 /d.o.f. = 22/23. The same fit without the additional term is shown by the dotted lines; the fit does not 
describe the data at low y/s well and x 2 /d.o.f. = 66/24 is found. One thus has evidence for additional contributions 
to this observables at low y/s probably behaving like \n(Q)/Q. 

Figure (right) shows results of simultaneous fits of as{mzo) and ao(2 GeV) to data for distributions of y® 3 at 
14 < y/s < 189 GeV. The solid lines show fits of the pure perturbative QCD prediction with an additional term of 
the form A 2 o/Q 2 inspired by Q. The fit describes the data at all cms energies and x 2 /d.o.f. = 71/106 is found. 
The dotted lines present the results of the same fits without the additional power correction term. The fit does not 
describe well the data at low cms energies, in particular at y/s = 14 and 22 GeV. One thus has evidence for additional 
terms probably behaving like 1/Q 2 at low y/s. At large cms energies the two types of fit are very similar. 
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Figure 4: Summary of results from power correction studies Q|. Data are from [9lllll.l2il.l23l. 



4. EVIDENCE FOR SOFT FREEZING 

One of the main assumptions of the DMW model of power corrections is, simply speaking, that the physical 
running strong coupling ag P remains finite around and below the Landau pole, where the perturbative evolution 
of the strong coupling according to the renormalisation group brakes down. With this assumption quantities like 
olq{im) = Jq 1 ag P (k)dk can be defined. 

Using hadronic decays of t leptons this question can be studied experimentally. The r lepton decays weakly via 
an intermediate (virtual) W + or W~ boson. When the W boson decays to quarks QCD effects have to be taken into 
account. The mass of the virtual W bosons m-yv < rn T — 1.777 GeV sets the scale for these effects which are rather 
large. The spectrum of invariant masses of hadronic final states from t lepton decays provides a direct experimental 
measurement of the scale given by the virtual m\y- 

Figure[S]left shows a measurement of the inclusive spectrum of invariant masses of non-strange hadronic final states 



s = ml by OPAL 26]. The quantity v(s) + a{s) = dR T /ds with R T (s) = T(r -» hj/ T )/F(r -> lv t v T ). The dashed 
line shows the "parton model prediction" for three colours. The solid line shows the perturbative QCD prediction 
for massless quarks. The line describes the correlated data points well; this implies that non-perturbative corrections 
are small for the inclusive spectrum. 

The analysis defines an effective charge a T (s) via R T (s) — i?°(l + a T (s)/7r). The relation between a r (s) and 
as(s) is known perturbatively. Figure (right) shows the variable a T (s) as a function of s extracted from the data 
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Figure 5: (left) Distributions of Bw measured at 14 < -*fs < 91 GeV as shown on the figure |2q|. Superimposed are fits 
perturbative QCD combined with power corrections with (solid and dashed lines) and without (dotted lines) additional terms as 
explained in the text, (right) Distributions of y^ 3 measured at 14 < y/s < 189 GeV as shown on the figure |2^|. Superimposed 
are fits of pure perturbative QCD (solid and dashed lines) or perturbative QCD combined with a power correction term as 
explained in the text. 
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Figure 6: (left) Spectrum of invariant masses of r lepton decays into non-strange hadronic final states |2g. The solid and the 
dashed line show predictions as indicated on the figure, (right) The effective charge a T determined from hadronic r decays as 
explained in the text as a function of s 17]. The lines show varios predictions for a T (s) as explained in the text. 



(shaded bands) and using different predictions (lines). The band from analysis of the inclusive spectrum labelled 
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"a T " stays fairly flat down to s ~ 1 GeV 2 and then rises steeply due to the p resonance. The other bands are from 
analysis of spectra for even or odd number of pions in the final state and have strong non-perturbative effects. The 
solid line labelled "/3 T 4 loop K4 = 25" corresponds to the most complete QCD prediction for the running of the 
effective charge including an estimate of 4th order term. This prediction describes the data well down to s ~ 1 GeV 2 
and stays finite at s = 0. One also observes that the NNLO prediction labelled "/3 r 3 loop" also describes the data 
at medium values of s and is finite at s = 0. ft is thus possible to consistently describe data for hadronic r decays 
down to s ~ 1 GeV 2 with an effective coupling which stays finite at s = 0. These observations are consistent with 
earlier work 0] • 



5. SUMMARY AND OUTLOOK 

In this review several topics were not covered due to a lack of results. One issue is the study of 4-jet event 
shape observables in the planar limit. The theoretical aspects and preliminary experimental results are discussed in 
contributions B001 and B002 to this conference. However, there are no published experimental studies which could 
be subject of this review, but hopefully there will be results on this important subject in the near future. Other 
models for the description of soft corrections to hard QCD processes such as shape functions, the single-dressed-gluon 
approximation (SDG) are subject of contributions to this conference but there are no detailed experimental studies. 
From an experimental point of view the DMW model is attractive, because it makes strong predictions about relations 
between observables and has only universal free parameters. Also, predictions actually exist for several observables 
such that direct comparisons are possible. 

The power corrections in the DMW model have been studied by several experimental groups as presented in this 
review. The picture is generally consistent with only a few notes: 

• The values for ao(2 GeV) from different observables are consistent within 20 to 30%. However, it has to be 
understood if this is a variation one expects from the theoretical uncertainties in the calculations. It is often 
claimed that the theoretical uncertainty on the Milan factor M of about 20% due to missing higher order terms 
could explain the scatter of ceo(2 GeV) results. However, since M is a universal factor it is not clear if it can 
explain differences between different observables. 

• The values for as (mz» ) from different observables are consistent with each other within the uncertainties of the 
perturbative QCD calculations as expected. However, direct comparison of the results from power correction 
analyses and analyses with hadronisation corrections derived from Monte Carlo simulation show a systematic 
bias for as(mzo) of about 6% (mean values) or 9% (distributions). This bias is comparable to or larger than 
the total uncertainties quoted for determinations of as( m z°) from event shape or jet observables in e + e~ 
annihilation. The bias is significantly larger than the hadronisation model uncertainties of less than 2% quoted 
for as (wz° ) determinations using LEP data Q] . 

In the future the studies of power corrections with c + c~ annihilation data would profit from further analyses using 
the JADE data in the intermediate energy region 14 < y/s < 44 GeV. In particular the EEC and 4-jet observables 
with or without requiring the planar limit would be attractive subjects. 

The impact of NNLO perturbative QCD calculations for 3-jet observables could also be significant. Firstly, the 
power correction calculations will have to adapted to the improved predictions. Secondly, with there is potential 
for high precision measurements of as(mz°) with reduced theoretical uncertainties. With selected observables with 
suppressed power corrections such as j/^ 3 one could obtain the matching reduced hadronisation uncertainties. 

Other approaches to understand soft corrections in terms of QCD should be studied in more detail. The connection 
with RS optimisation discussed in contributions T001 and T002 is also an interesting new area of research. 
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